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Abstract 

We show that the two dimensional Ising model is complete, in the sense that the partition func- 
tion of any lattice model on any graph is equal to the partition function of the 2D Ising model with 
complex coupling. The latter model has all its spin-spin coupling equal to i? and all the parame- 
ters of the original model are contained in the local magnetic fields of the Ising model. This result 
has already been derived by using techniques from quantum information theory and by exploiting 
the universality of cluster states. Here we do not use the quantum formalism and hence make the 
completeness result accessible to a wide audience. Furthermore our method has the advantage of 
being algorithmic in nature so that by following a set of simple graphical transformations, one is 
able to transform any discrete lattice model to an Ising model defined on a (poly normally) larger 
2D lattice. 

1 Introduction 

Despite the fact that foundations of statistical mechanics were laid down in the nineteen century, the 
understanding that a single partition function can describe different phases of matter, is rather recent 
and indeed as late as 1930's, there was not a consensus among physicists that one single partition 
function can describe ice, water and vapor (TJ or can give a sharp phase transition which we observe 
in ferromagnetic materials. The works of Kramers, Wannier, Onsager lEJEHU and others gradually 
established beyond doubt that in the thermodynamic limit, singular behavior and phase transition can 
arise from a single partition function based on a single model Hamiltonian. For example the Ising 
Hamiltonian can describe both the ordered phase of a ferromagnet and the disordered phase of a 
paramagnet. However, the study of second order phase transition over the past decades, revealed an 
interesting feature which is called universality. It means that the near the point of second order phase 
transition, the details of the model Hamiltonian do not matter and only some general properties like 
symmetries and dimension of the order parameter are important. 
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Our current understanding of this important phenomena, that is the independence of the most 
salient features of phase transition from the microscopic details is based on Renormalization Group 
(RG) and effective theories. Under renormalization group (change of scale), different microscopic 
theories tend toward an effective theory which describes the long range characteristics of a whole 
class of models, irrespective of their microscopic details. Therefore if we imagine an infinite dimen- 
sional space in which all forms of Hamiltonians describing the interactions of Ising type variables 
are defined, there are some special points in this space, to which other models tend toward the RG 
flow. These fixed points essentially capture the most important properties of phase transitions of all 
the models along the RG trajectories. 

Quite recently a new notation, different from universality, has emerged in the study of classical 
statistical mechanical models which is appropriately called completeness J6] |7] |9l \W\ HU [T2) . The 
idea of completeness asserts that certain statistical mechanical models are complete in the sense that 
if the partition function of these complete models are solved, then the partition function of any other 
model can be also solved with a polynomial computational overhead. This new insight certainly is 
different from the notion of universality which applies only to the long distance limit of different the- 
ories. It takes into account the microscopic details of different models and yet claims that all these 
models are inherited from one single model. It is reminiscent of the concept of complete problems in 
computability theory 1251 26). Therefore, in the same space of coupling constants mentioned above, 
there are certain special points, from which a whole class of other models can be derived, in the sense 
that the partition function of the latter is exactly equal to the partition function of the complete model 
if we formulate the complete model on a sufficiently large lattice. 

Up until now many models have been shown to be complete J6] |7] [H |9] [TO) . In this paper we are 
concerned with the Ising model which was the first to be proved complete |6|. It is defined on a 2D 
rectangular lattice, with the following Hamiltonian 

where Si = ±1, the hi's are local magnetic fields which can be complex and all the interactions are 
nearest neighbors. 

It is true that the model has no physical realization due to the complex fields and couplings, how- 
ever it is quite significant that any other model, regardless of the nature of its statistical variables, 
dimensionality and even the geometry of its lattice, and the nature of its interactions is just an instant 
of this general computational problem. All the details of the original model, the couplings, the mag- 
netic fields, the connectivity and geometry of the lattice are contained in the size of the 2D lattice and 
the inhomogeneous complex magnetic fields. 

Almost all of these completeness results J6] H HJ |9] [10] QT| [T2) were derived by merging of ideas 
from statistical mechanics and quantum information theory lTT4l [T31 fT6l rT71 [181 [T9l l20l ETl l22l l23l 
24]. The new paradigm of Measurement-based Quantum Computation (MQC) and the universality of 
cluster states for MQC j27]|28l|29l|30l|3l][32[33[Sl[32 essential in 

these derivations. 

The starting point was the observation in J5] that the partition function of any given discrete model 
can be written as a scalar product, 

Z G (J) = (a|* G ), (2) 

where (a| is a product state encoding all the coupling constants J and \^g) is an entangled graph 
state, defined on the vertices and edges of a graph, and encoding the geometry of the lattice. The core 
concept of completeness is the fact that the 2D cluster state is universal. More concretely we know 
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that the graph state |$g) corresponding to a graph G can be obtained from an appropriate cluster 
state l^n) corresponding to a rectangular lattice (denoted by □), through a set of adaptive single- 
qubit measurements Ai. The measurements being single-qubit can be formally written as \<xm) (&m\, 
where |q;m) is a product state encoding the bases and the results which have been measured. Thus the 
totality of measurements Ai transforms the cluster state as follows 

!*□) — -> \a M )(a M \y n ) = \a M ) ® |* G ), (3) 

which shows that after disregarding the states of the measured qubits \(xm), what is left is an appro- 
priate graph state 

|* G ) - (a M |*n>. (4) 

Combination of this relation with (O, leads to the completeness result mentioned above. That is one 
writes 

Z G {J) = (a\9 G ) = (a, a M \9n), (5) 

and notes that (a, a^/| now encodes a set of generally inhomogeneous pattern of interactions on the 
cluster state. Therefore one has 

Z G (J) = Z n (J,J'). (6) 

However, there are two shortcomings in the above derivation. First it is an existence rather than a 
constructive proof, and one does not learn how to reduce any given statistical model to the 2D Ising 
model, at least not in a simple way which we intend to introduce. Second, it is not accessible to those 
who are not familiar with the concepts and techniques of quantum information theory, notably that of 
measurement-based quantum computation. While the result is a general one in statistical mechanics 
with potential intersect across many disciplines, its derivation is somewhat specialized. 

It is therefore the aim of our paper to provide a proof for this completeness result which is both 
constructive and accessible to a large audience and readership. We show how any model can be re- 
duced to an instance of 2D Ising model through a set of simple graphical transformation steps. As a 
by-product we also re-derive the well-known duality relations between various Ising models, in a way 
which is different from the high and low temperature expansions. Finally our proof does not require 
any basic knowledge of quantum information theory or MQC or even quantum mechanics. We hope 
that by providing this proof, we will make this important result accessible to a wide audience who can 
further explore the far-reaching consequence of the completeness of 2D Ising model. 

The structure of this paper is as follows: In section 10, we first review the basic result that any 
discrete statistical mechanical model can be expressed in terms of Ising variables. We then show how 
the interaction terms of any such model can be transformed to the local magnetic fields of an Ising 
model whose interaction strength is fixed at the value of i\. In section (O we apply this idea to a few 
well-known examples and use it to derive the well-known results on duality in a way, which is inde- 
pendent from the high and low temperature expansions. In section we come to the main result of 
our paper, which is the proof of completeness of 2D Ising model. We introduce basic transformation 
rules and express them in terms of diagrams (with further examples in an appendix) which allows us to 
transform any Ising model on any kind of graph (planar or not), to the 2D Ising model on rectangular 
lattices. The paper ends with a discussion and the above mentioned appendix. 



2 Transformation in lattice models 

Consider an arbitrary discrete lattice model, like Ising, Potts, vertex or face model |[T3l . The discrete 
variables of this lattice model need not be bi-valued like the Ising model, and the lattice need not be 
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regular. It has been shown in (7) that the Hamiltonian of such a model can always be re-written in 
terms of two-valued spin variables , in the form 

H = ^ JiSi + ^ JijStSj + ^2 J'i]kSiS S k + (7) 

where Si = {1, —1} is a spin variable and J i5 Jy-, .. are coupling constants. Therefore the price that 
one pays for converting any type of variable to the Ising variable is to admit few-body interactions 
in the Hamiltonian in addition to two-body interactions. As two simple examples, consider the the 
three-state and four- state Potts model, whose variables q take 3 or 4 different variables respectively, 
and their Hamiltonians are given by: 

ff=-^Ew (8) 

(id) 

As shown in appendix A, by suitable encodings of the Potts variables into a a pair of Ising variables 
(Si, S'i), we can turn these models into Ising models on lattices of 2N sites. 

These two examples, illustrate the general idea of ^} that by encoding any kind of discrete vari- 
ables into a number of Ising variables, one can turn any statistical model into an Ising model, albeit 
with complicated multi-particle interactions. 

Having this result, it is now enough to consider Ising model on arbitrary graphs with arbitrary 
patterns of interactions. The Hamiltonian of such a model is written as 

H = ^ JiSi + ^ JijSiSj + ^ Jijk.SiSjSk + ■ ■ ■ , (9) 

where Jj, Jy and Jy*. are the strengths of one, two and three body interactions and • • • stand for 
interactions of higher number of spins. We can now turn this Hamiltonian into an equivalent one, 
containing only on-site interactions or local magnetic fields. The procedure is to define new spin 
variables as 

Sij = SiSj , Sijk = SiSjSk, ■ ■ ■ 5 (10) 
which turn the above Hamiltonian into 

// X JA + Yl J v s » + S J mSijh + ■■■ ■ (ID 

Note that the variables Sij are assigned to links, the variables Sijk to triangles and other variables to 
higher simplexes of the lattice. 

If all these new variables were independent, then the new Hamiltonian would have been that of 
a collection of spins S a € {Si, Sij, Sijk, • • • } in local magnetic fields h a € {Ji, Jij, Jijk} with a 
trivial partition function. However, the price that we should pay for this change of variable is that 
the new variables are not independent. In fact as we see, this price has the extra merit of paving the 
way for the proof of completeness result. There are many constraints between the new variables, for 
example we have 

SiSijSj = 1 , SiSjkSijk = 1, • • • . (12) 

In order to calculate the partition function in a correct way, one should take into account all these 
constraints. Therefore the general form of the partition function will be 

Z = Ar £ e -^^Tp ( -Qs Q ,l). (13) 

{S a } C a£C 
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Figure 1: (Color Online) Graphical representation of lemma ( fT4b . A white circle denotes a spin which 
is summed over. The values on spins indicate the local magnetic fields. Any link indicates a coupling 
with value m is the number of black spins. 

In the above formula, C denote a constraint and a G C stands for the spin variables taking part in a 
constraint C, and A is a numerical factor which takes of overcounting. 

The basic idea is now to rewrite the constraints, which contain product of spins, in such a way to 
turn them into a form containing only two-body interactions, no matter how complicated the constraint 
is. To do this we introduce the following lemma: 

Lemma : For each constraint, we can add one single spin variable Sq and write it as 

S(S 1 S 2 ---S m = l) = l £ e f(i-s»)(™-a-&-i). ( i4) 

So=l,-l 

When inserted into the partition function, this produces extra magnetic fields —i\ on each spins Si, a 
magnetic field of — m^f on the extra spin Sq and two-body interactions of strength ^ between pairs 
of spins. The proof of lemma is very simple: 

Proof : We use the variables a <G {0, 1} defined as S = (— 1) CT = e J7rcr to write the constraint as 

6(S 1 S 2 ---S m ,l)=8(a 1 + a 2 + a n =0) = ± £ e *™» («+«+•••+*»). (15) 

" er o =0,l 

Now going to the original Ising variables through the relation a = converts this last relation to 
equation ( [Pil l, hence the lemma is proved. Lemma dT4b is depicted graphically in figure ([T). For this 
and the other identities which are depicted in other figures, we will use a convention which we gather 
as follows: 

Notations and Conventions: In all the equations and figures in this section, the Boltzman factor 
— ft = ^--L has been absorbed in the Hamiltonian, so H means —f}H in all the equations. The param- 
eter 7 has the fixed value of if. A spin which is summed over, is depicted by a white circle. A gray 
circle means a spin whose value has been fixed to S — 1. Every link indicates an interaction with a 
coupling strength of ij between its endpoints. The value written near any white circle is the value of 
the local magnetic field on this spin. When a + or — sign comes on the right hand side of a numerical 
value on spin, it means that the effect of summation over the white spin is, the addition or subtraction 
of that value to the local magnetic field of that spin. 

Using this identity, the partition function of any model will be equivalent to the partition function 
of an Ising model on a complex graph in which the coupling constants of all bilinear couplings is the 
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Figure 2: (Color Online) An Ising model on a triangular lattice. On each link of the lattice, we put a 
new spin Sij := SiSj. These new spins are constrained, in that the product of all of them around a 
triangle is equal to 1 . This is accounted for by inserting a white spin inside each triangle and using 
lemma (TT~4-b . The result is a hexagonal bi-partite lattice. Summing over the black spins will lead to 
duality between Ising models on triangular and hexagonal lattices. See subsection (13.4b . 

constant value of ij. It is clear that for any model which contain interactions of 2-body, 3-body up 
to a finite k — body interactions, the number of constraints is polynomially dependent on the number 
of vertices of the original graph and hence the size of the final Ising graph is only polynomially larger 
than the original graph. Moreover our construction shows a specific algorithm for turning the original 
discrete model to an Ising model with only ij interactions. We still have to show that the new Ising 
model on the complex graph which may well be non-planar, is equivalent to a 2D Ising model. We will 
prove this completeness result in section (0]). However before doing this, it is instructive to present a 
few examples to show how the above transformation maps some well known models to other ones. 

3 Examples 

In this section we consider a few examples to show how the transformation (fT~4T >. makes a discrete 
model equivalent to an Ising model on a different lattice. The first two examples deal with Ising 
model with two-body interactions on triangular and hexagonal models respectively, while the third 
one deals with a new model introduced recently in l44ll45l with three body interactions. 



3.1 Ising model on triangular lattice 

Consider an Ising model without magnetic field on a triangular lattice with free boundary conditions. 
The partition function of this model is 



On each edge e := (i, j) of the lattice, which connects the vertices i and j, we insert a new variables 
S e = SiSj . These new variables are not independent. In fact for each triangle of the lattice with 
edges e\ = e2 = (j, k) and = (k, i), we have S ei S e2 S e3 = 1. Implementing this constraint 

via the relation ( TPfl i or the corresponding figure (Q]), we arrive at a bi-partite hexagonal lattice, shown 
in figure (0, with two kind of vertices, which we designate as black (B) and white (W). The black 
circles designate the spins S e put on the edges of the original lattice, and the white circles designate 




(16) 
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the new spins which take care of the constraint through equation (TPfl i. From figure 0, we find that 
Hamiltonian of the new model is as follows: 

H hex = {j-^)Y,Si-^ E £ SiS - (17) 

ieB j£WB i£B,j£W 

Taking into account the multiplicative factor for each triangle in figure 0, we arrive at the following 
relation between partition functions: 

Z A (J) = (^e 3i ^Z hex (J), (18) 
where | A| is the number of triangles in the original lattice. 

3.2 Ising model on hexagonal lattice 

Consider now the Ising model on a hexagonal lattice as shown in figure (0. The Hamiltonian is 

H hex = J J2 S i S h (19) 

itihex 

where hex means all the vertices of the hexagonal lattice. Again we do the same transformation as 
before, that is on any edge e = connecting the two vertices i and j which allocate the spin 

variables Si and Sj, we put a new spin variable S e := SiSj. The original coupling constant J in 
the Hamiltonian (fT9l l will now play the role of a magnetic field in the new lattice. However the new 
variables are not independent and indeed on for any hexagon in the lattice, the product of all the spin 
variables around the edges is equal to 1 . Hence using (fl4l i we find that the model is equivalent to a 
new triangular lattice, as shown in figure 0. Again the new lattice is bi-partite, with two sets of sites 
denoted as black (B) and white (W) as in the previous example. The Hamiltonian of the new model 
is: 

F A = (J-i|)^5 i -^E^'+ i I £ Mi- (20) 

ies jew ieBjew 

In view of the relation in figure (0, the relation of the partition functions will be 

Z hex (J) = (\e 3l %)^Z A (J), (21) 
where \hex\ is the number of hexagons in the original lattice. 



3.3 An Ising model with three-body interaction 

The final example is a kind of Ising model on triangular lattice with three body interactions. This 
example was studied in [43] in the context of topological color codes, it is in a different universality 
class than the usual Ising universality class 03) ■ The Hamiltonian of this model is 



H 3 A = J SiSjSk, (22) 
(i,j,k) 

where by fc), we mean the vertices belonging to a single triangle. In this case, the new spin 
variables are assigned to each triangle (plaquette), as St ■— SiSjSk, where i, j, and k are the vertices 
of the triangle t. The coupling constant J plays the role of local magnetic field on each plaquette. 
The new spin variables satisfy the constraint that ]X St = 1, where by Y[' we mean the six triangles 
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Figure 3: (Color Online) An Ising model on a hexagonal lattice. On each link of the lattice, we put a 
new spin Sij := SiSj. These new spins are constrained, in that the product of all of them around a 
hexagon is equal to 1 . This is accounted for by inserting a white spin inside each hexagon and using 
lemma (TPfl) . The result is a triangular bi-partite lattice. Summing over the black spins will lead to 
duality between Ising models on hexagonal and triangular Ising lattices. See subsection (13.4b - 

sharing one single vertex in the center. The new lattice is shown in figure (01, and is seen to be 
triangular again, which we denote by A*, with the following Hamiltonian: 

ffi. = (J-i^)$>-3i££s,- +i I J2 S * S i- ( 23 ) 

i£B jeW ieBJeW 

Taking into account the multiplicative factors in figure we find the following relations between 
the partition functions: 

Z A (J) = (\e 3 ^)^Z A ,(J), (24) 

where \t\ is the number of triangles in the original lattice. Note that the new lattice has only two-body 
interactions. 



3.4 Duality 

As the reader may have noticed, the graphs obtained for the above examples, are all bi-partite graphs. 
This allows us to do a partial summation over the subgraphs, comprised of black spins and obtain an 
equivalent Ising model. This then leads to the well-known duality relations. As it is evident from 
figure ([3]), doing such a partial summation on the final Ising graphs which were derived from the 
hexagonal and triangular Ising models, leaves us with Ising models on triangular and hexagonal lat- 
tices respectively, but with different couplings. To explicitly demonstrate this duality, we detail it only 
for the rectangular lattice which is known to be self-dual. Repeating the analysis for other kinds of 
models is straightforward. 

Figure © shows the 2D Ising model. The light-blue circles are the spins on the original lattice, 
interacting with the couplings J. The black spins are the new spins which are assigned to each link 
as in ( [Tol l and the constraint between these new spins around each plaquette, is taken care of by 
addition of a single white spin in each plaquette. The new bipartite lattice is governed by the following 
Hamiltonian: 

f}H = (J- i ^) y £s i + i* , ES j + i ^ ]T SiSj, (25) 
ieB jew ieB,je{weNi] 



Figure 4: (Color Online) An Ising model on a triangular lattice with three-body interactions. On each 
face of the lattice, we put a new spin Sijk '■= SiSjSk- These new spins are constrained, in that the 
product of all of them around a hexagon is equal to 1 . This is accounted for by inserting a white spin 
inside each hexagon and using lemma (fT4l . The result is a triangular bi-partite lattice. 



where the following relation exists between the original partition function and the new one: 

Z2DIsing{J) = % 2D I sing ( 26 ) 

Here N is the number of spins of the rectangular lattice and Z* is the partition function of the equiv- 
alent model. If we now sum over the black spins, we obtain a new Ising model comprised entirely of 
white spins. To find the Hamiltonian of the new spin model, a simple calculation is in order. 

Consider a black spin So and it's two neighboring white spins, Si and S2 ■ Summing this part of 
the partition function we find 

j2 e ^s (s 1+ s 2 )+(J-^)So = 2cosh(^(5i +S 2 ) + (J- y)). (27) 

So 

After a simple algebra, the right hand side can be rewritten in the form 

cosh^Si + S 2 ) + (J - y )) = Ae KSlS *+ hlSl+h » s > (28) 
where the parameters are required to satisfy the following relations: 

Ae K+h 1 + h 2 = cogh (J) ; 

Ae -K+ hl -h 2 = _j sinh (J) ; 
Ae ~K~h 1+ h 2 = _ isinil (J-) } 
Ae K-hi-hz = _ cosh (J), 



the solution of which is 



hi = h 2 = — 
A 2 = -isinh(2J) 

e~ 2K = tanh{J). (29) 
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Figure 5: (Color Online) The Ising model on a rectangular lattice and its dual. The light-blue circles 
are the spins on the original lattice, interacting with the couplings J. The black spins are the new 
spins which are assigned to each link as in (fTOb and the constraint between these new spins around 
each plaquette, is taken care of by addition of a single white spin in each plaquette. 



The last relation is the well-known duality relation in the 2D Ising model, which is usually obtained 
by comparing the high and low temperature expansion of the Ising model partition function. Here it 
is derived without any reference to high or low temperature expansions. Putting all this together, we 
find 

Z m ( J ) = 2N sinh(2J) Ar Z 2DIsing (X), (30) 

where N is the number of spins of the 2D lattice. This is exactly the duality relation of the 2D Ising 
model. This type of analysis can be repeated for triangular and hexagonal lattice models with ensuing 
duality relations between them. These examples show how the mapping works and how every type 
of interaction in the original Hamiltonian turns into a magnetic field and how the constraints turn 
into complex two-body interactions of ij in the new model. We have deliberately considered simple 
examples in which the constraint do not produce complex non-planar graphs. In general this can 
happen and the new graph can be topologically complex. Nevertheless we show that by a series of 
simple transformations, the new model can be turned into a 2D Ising model. Therefore we provide 
an algorithmic approach for transforming any discrete model to the 2D Ising model, and in this way 
we provide a different proof of completeness of 2D Ising model. Compared to the existence proof (6j 
based on universality of cluster states, this new proof is constructive in nature. 



4 Completeness of the Ising model on 2D square lattices 

In previous sections we have shown that any discrete statistical mechanical model is equivalent to 
an Ising model with complex two-body interactions of ^ where the effect of the original coupling 
constants is transferred to the local magnetic fields on the spins. Depending on the original statistical 
mechanical model, the graph of the equivalent Ising model may be complicated and usually it may be 
a non-planar graph. We now introduce the steps by which one can transform such an Ising model to 
one on a polynomially larger two-dimensional rectangular lattice. The steps or the rules go with the 
name of uniformizing and flattening steps respectively. We describe them in that order. The notations 
and conventions introduced before, apply to all which follows. 
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Figure 6: (Color Online) Identity OTb : the white spin has zero magnetic field and has been summed 
over. This allows us to increase the degree of any vertex as required. Note that the magnetic field of 
the black spin does not change, hence the notation 0+ on this spin. 




Figure 7: (Color Online) Identity d32l i. The white spin is in magnetic field «'| and is summed over. 
The result is merging of the two neighboring spins. Reading the diagram from left to right, it shows 
how a high degree can be split to two vertices of lover degree. 



4.1 Uniformizing the Graph 

As shown in section (f2]i and figure (Q~|i, for any constraint between a set of spins {Si, S2, ■ • • S m } on 
vertices {1, 2, • • • m}, we have to add a new vertex with spin variable S, which links connecting this 
new vertex to all the vertices 1 through m. This may produce a complex graph with vertices with 
various degrees of connectivity. The uniformizing rule allows us to transform the graph so that the 
degree of all its vertices, modulo those on the boundary are equal to four. We do this in a few steps. 
First, we note the following simple identity: 

This identity depicted in figure ((6]), shows that we can add any isolated vertex with ani| interaction 
to an existing vertex without any effect except a multiplicative factor of \J~2 which can be taken care 
of in the final partition function. 

To reduce the degree of vertices down to four, we present the following 

Lemma: For Ising variables (taking values ±1), the following identity holds: 

S SiS2 = l£y§So+iiS (s 1+ s 2 )_ (32) 

To verify it, we note that the right hand side is equal to 

_ I ( ie «? (Si+&) _ ie -ii = shi l( Sl + S2) = dsiS2 _ (33) 

This identity depicted graphically in figure © allows us to reduce the degree of any vertex down to the 
value four. Actually as this figure shows, the degree can be reduced to three, however we need degree 
four for the flattening step which will be explained later. Finally we note that we can insert vertices 
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Figure 8: (Color Online) Identity d34b . Insertion of one vertex on a link: The white spin with a 
magnetic field ij can be inserted on a link, reducing the magnetic fields on the endpoints by i^. 
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Figure 9: (Color Online) Repeating the basic insertion in figure ([8]). a) two vertices and (b) three 
vertices are inserted on the link connecting the black vertices. 

into already existing links and plaquettes. The insertion in a link follows the following identity: 

So+ 1 f L So(5 , i+S2) _ _|_ j\g i j L (Si+S2+SiS2) (34) 



So 



This identity is diagrammatically shown in figure (8). It shows that a spin with magnetic field ij can 
always be inserted in a link with the effect of reducing the magnetic fields on the end vertices of the 
link by ij. Repeating this insertion once and twice we arrive at figure (|9}a) and (|9]-b) respectively, 
with a similar pattern for larger number of insertions. 

Insertion into a plaquette is achieved simply by fixing the value of that spin to 1 followed by 
reducing the magnetic fields on all the spins on the vertices of the plaquette by ij. This is shown in 
figure ( [Tol l and the verification of the identity of the two sides of the figure is trivial. 

Therefore with these identities we can uniformize all the degrees of the graph and turn it into a 
polynomially larger graph with degree four for all vertices. Some examples for working with these 
steps are given in the appendix B. 

4.2 Flattening the Graph 

The graph obtained by implementing the constraints, may well become a non-planar graph. In this 
subsection we show that any such graph can be made a planar graph by adding a number of vertices 
of the same order as the number of vertices in the original graph. Again we present a useful identity 
in the form of the following 

Lemma: Let gq be a qubit variable, taking values and 1. Then for any integer x, the following 
identity holds: 

^- eif (ao+^) 2 = ( 1+ { y (35) 
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Figure 10: (Color Online) Insertion of a spin into a plaquette. The value of the central spin is fixed 
to 1 and the magnetic fields of all the spins connected to it are reduced by if, with no change in the 
partition function. 

The proof is by simple substitution and by noting that the square of any two consecutive integers 
are of the form 4fc and 4fc + 1. From this we can then prove the following: 

Lemma: Let <jq and Oi be qubit variables. Then the following identity holds: 

(70=0,1 

where K is a subset of the edges which are the neighbors of erg and K' is the complement of that 
set. To prove this we complete the square in d36l l and use lemma (l35l l with x = YllLi <J i- More 
specifically we write the left hand side of (l36l l as 

LHSof(0ni)= e^^+^e-^^^-™^-.^*" 7 ^. (37) 

<T =0,1 

Using lemma d35l l and noting that e l1Ty = e~ mY for any integer, we arrive at ( l36l >. 

The lemma ( T36T ) which is reminiscent of what is called local complementation rule in the context 
of measurement-based quantum computation, plays a basic role in flattening every diagram. First note 
that while the form of the lemma is general in terms of qubit variables (a = 0, 1), when expressed 
in terms of Ising variable (S = 1, —1), it may not have a simple form, independent of the details of 
the connections of the graph. However we need this lemma only for a special simple case. This is the 
graph of figure ( TTTb . For this simple case the substitution of qubit variables to Ising variables, gives 
the following simple identity 

^ e -jf5 , o+ifS , o(Si+S2+S 3 +S4)+if(Si+S 3 )(S2+S4) _ ^ _ ^gif (Si+S 2 +S 3 +S 4 )+if {S 1 S 3 +S 2 S i ) 
Sa 

(38) 

This identity is shown in graphical form in figure (fTTT i and it allows us to remove all the over-crossings 
and hence flatten all the diagrams. 

Therefore by a combination of insertion into the links and faces and removing all the crossings, 
we can turn the arbitrary graph which is obtained after imposing the constraint into an Ising model on 
a 2D rectangular lattice. The number of vertices which should be added grows polynomially with the 
number of spins in the original statistical model, provided that the range of interactions is bounded by 
some integer independent of the size of the lattice. In this way it is proved that any discrete statistical 
model is equivalent to a special kind of rectangular 2D Ising model, with the following Hamiltonian: 

H = J2hiSi + ijJ2 SiS ^ (39) 

1 
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Figure 11: (Color Online) Identity d38b . The white spin is in magnetic field — ij and is summed over. 
The result is complementing the interaction links, as shown, with the addition of magnetic fields if 
to the fields already present in the neighbors. This diagrams allows us to remove all the crossings and 
flatten the arbitrary graph. 

that is all the neighboring spins interact with each other by the coupling constant and all the 
original coupling constant and the geometry of the original lattice is reflected into the inhomogeneous 
magnetic fields on the local spins. 

5 Discussion 

The Ising model in two dimensional rectangular lattices, i.e. 2D Ising model, is a complete model 
in the sense that any the partition function of any other statistical mechanical model with discrete 
type variables, is equivalent to a 2D Ising model on a sufficiently large rectangular lattice. This Ising 
model, whose Hamiltonian is in the form of (13 91) . is special in the sense that all its coupling constants 
are equal to the complex value of ij and all the coupling constants of the original model and all the 
details of the geometry of the original lattice, and the pattern of interactions are reflected into the local 
magnetic fields of the Ising model. This completeness, was first proved J6) by using concepts and 
techniques from quantum information theory and was based on the universality of cluster states. We 
have now proved this important result, by using independent and general concepts and methods which 
are accessible to a wide audience of researchers across many disciplines. In fact our method, not 
only unravels the essence of this completeness from a new perspective, it also provides an algorithmic 
approach through definite and simple graphical steps for reducing any other model to an Ising model 
of the above type. The universal form of the Ising model as shown in Eq. d39l ), may have far reaching 
consequences. It shows that all the various models in different universality classes are different only 
in the number of sites and the pattern of local magnetic fields on this 2D lattice, and not on the 
nature of discrete variables, nor on the pattern of interactions. In particular the well-known criteria of 
dimensionality of the system and symmetry of the order parameter have to be re-interpreted in terms 
of the complete 2D Ising model. 
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6 Appendix A 

In this appendix we provide two simple examples to illustrate the assertion in section (O that the 
discrete variables of any kind of statistical model can always be expressed in terms of Ising variables. 
The examples are the four-state and three-state Potts models. 
Consider 1 dimensional four-state Potts model with the following Hamiltonian: 

(hi) 
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Figure 13: (Color Online) Example 1: Insertion of a vertex on a link, as in (|S), turns a triangle to a 
rectangle, with appropriate changes in the magnetic fields of the adjacent vertices. 

where each qi takes four different values {0,1,2,3}. We can encode a variable qi into two Ising 
variables (Si, 5|) and hence turn the one-dimensional chain into a ladder as in figure (fT2l and use the 
encoding 

= (0,0), 1 = (0,1), 2 = (1,0), 3= (1,1). 
A simple calculation then shows that 

K,n = |(! + S iSj + S 'iS'i + SiSjS'iS^). (41) 
For the three-state Potts model, we use the following encoding 

= (0,0), 1 = (0,1) -(1,0), 2 = (1,1), 

where only each Potts value is mapped to the equivalence class of the Ising pair, under symmetry. In 
this it can be verified that 

3 1 3 

3qi,qj = g ~t~ g( SiS^ ^j^j "F SiSj -F S^Sj + SiSj + S^Sj) -F —SiS^SjSj. (42) 

Therefore both models can be expressed in terms of Ising variables. 

7 Appendix B 

In this short appendix we show how the uniformizing rule is to be used via a few simple examples. 
We present these examples only to show the working of the basic identities. Clearly converting a very 
complex graph requires many steps to be taken, but the nature of them all is the same as a kind of 
uniformization and flattening. We consider three simple examples: 

Example 1: 

Consider a simple triangle as shown in figure (|F3T >. Insertion into a link, turns this into a rectangular 
plaquette, the simplest 2D Ising model. 

Example 2: 

Consider a pentagon as shown in figure ( TPfl i. Three link-insertions followed by a face-insertion turns 
this into the rectangular lattice on the right hand side. 

Example 3: 

Consider a crossing as in figure (TT5l ). Flattening this according to (TTTb . followed by insertions on the 
four links turns this into a rectangular lattice with four squares. 
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Figure 14: (Color Online) Example 2: Insertion of three vertices on the links of a pentagon, as in 
©, followed by an insertion of a vertex into the plaquette as in (TTOb , turns a pentagon into a small 
rectangular lattice. 




Figure 15: (Color Online) Example 3: By local complementation and insertion of four vertices on the 
links, we can remove a crossing. 
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